Melting of the vortex lattice in high T c superconductors 
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The precise measurements of vortex melting point towards 
a need to develop a quantitative theoretical description of 
thermal fluctuations in vortex matter. To tackle the difficult 
problem of melting, the description of both the solid and the 
liquid phase should reach the precision level well below 1%. 
Such a theory in the framework of the Ginzburg - Landau 
approach is presented. The melting line location is deter- 
mined and magnetization, specific heat jumps along it are 
calculated. We find that the magnetization in the liquid is 
larger than that in the solid by f .8% regardless of the melt- 
ing temperature, while the specific heat jump is about 6% 
and slowly decreases with temperature. The magnetization 
curves agree with experimental results on YBCO and Monte 
Carlo simulations. 



A magnetic field generates an array of vortices in type 
II superconductors. The vortices strongly interact with 
each other forming highly correlated configurations such 
as the vortex lattice. In high T c cuprates at relatively 
high temperatures, vortices move and vibrate due to ther- 
mal fluctuations to the extent that the lattice can melt 
becoming a " vortex liquid" [1]. Several recent remarkable 
experiments clearly determined that the vortex lattice 
melting in high T c superconductors is a first order phase 
transition with magnetization jumps [2] and spikes in 
specific heat [3]. The magnetization and entropy jumps 
were measured using diverse techniques such as local Hall 
probes, SQUID [4,5], torque magnetometry [6,7] and in- 
tegration of the specific heat spike [3] . Related investiga- 
tions indicate that, in addition to the spike, there is also a 
jump in specific heat [3,8]. In spite of those precise mea- 
surements of vortex melting, a quantitative theoretical 
description of vortex melting is still lacking. We present 
such a theory in the framework of the Ginzburg - Landau 
approach in this work. 

Thermal fluctuations in vortex matter have attracted 
great attention since the high Tc superconductors were 
discovered over a decade ago. In highly anisotropic su- 
perconductors like BSCCO near the melting point vor- 
tices are quite well separated and the system can be ap- 
proximated by an array of pointlike objects [1]. In less 
anisotropic ones like YBCO, near the melting point the 
vortices overlap and one has to use the Ginzburg - Lan- 



dau (GL) model. The model is defined by the following 
free energy: 
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where A = {By. 0) describes magnetic field (consid- 
ered constant and nonfluctuating, see below) in Lan- 
dau gauge and covariant derivative is defined by D = 
V - «i^A,$ = (e* = 2e).Whcn k = A/£ is large 
(greater than 10 ) where A magnetic penetration depth 
and £ coherence length, the magnetic field can be consid- 
ered constant and nonfluctuating (an excellent approxi- 
mation in the region studied). Statistical physics is de- 
scribed by the statistical sum 

It accurately describes thermal fluctuations in the range 
of relatively large magnetic fields (H >> H c i = 100G 
in YBCO) and temperatures near T c (70K to llCK" 
in YBCO ). Near T c , a(T) can be approximated bt 
aT c (l — t). In this paper, we will consider only high k 
and temperatures near T c . The model is, however, highly 
nontrivial even within the lowest Landau level (LLL) ap- 
proximation. In this approximation, which is valid when 
the magnetic field is high, only lowest Landau level mode 
is retained and the free energy simplifies (after rescaling) 
[9]: 
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The simplified model has only one parameter, the di- 

( b \~ 2/3 

mensionless scaled temperature 
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The dimensionless Ginzburg number Gi characterizing 
the importance of thermal fluctuations is \ " 

and the anisotropy parameter 7 is \fm c /m a i,. 

The LLL GL model was studied by a variety of differ- 
ent nonperturbative analytical methods. Among them 
the density functional [10], l/N [11], elasticity theory 
[12] and others [13]. The model was also studied numeri- 
cally in both 3D [14] and 2D [15,16] However we will not 
discuss those approaches in this paper. 

While applying the renormalization group (RG) on the 
one loop level to this model, Brezin, Nelson and Thiaville 
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[17] found no fixed points of the (functional) RG equa- 
tions and thus concluded that the transition to the solid 
phase is not continuous. The RG method therefore can- 
not provide a quantitative theory of the melting transi- 
tion. Two perturbative approaches were developed and 
greatly improved recently to describe the solid phase and 
liquid phases, respectively. The perturbative approach 
on the liquid side was pioneered by Ruggeri and Thou- 
less [18], who developed an expansion in which all the 
"bubble" diagrams are resummed. Unfortunately they 
found that the series are asymptotic and, although first 
few terms provide accurate results at very high tempera- 
tures, the series become inapplicable for ax less than —2 
which is quite far above the melting line (located around 
ar = —10). We obtained recently an optimized gaus- 
sian series [9] which is convergent rather than asymptotic 
with radius of convergence of ar ~ —5, but still above 
the melting line. 

On the solid side, Eilenberger [19] calculated the fluc- 
tuations spectrum around Abrikosov's mean field solu- 
tion. Maki and Takayama [20] noticed that the vortex 
lattice phonon modes are softer than those of the usual 
acoustic phonons in atomic crystals and this leads to in- 
frared (IR) divergences in certain quantities. This was 
initially interpreted as the destruction of the vortex solid 
by thermal fluctuations and the perturbation theory was 
abandoned. However the divergences resemble the "spu- 
rious" IR divergences in the critical phenomena theory. 
A recent analysis demonstrated that all these IR diver- 
gences cancel in physical quantities [21] and the series 
therefore are reliable. The two loop calculation was per- 
formed, so that the LLL GL theory on the solid side is 
now precise enough even for the description of melting. 

However, on the liquid side, a theory for cit ^ —5 is 
required. Moreover this theory should be extremely pre- 
cise since the internal energies of the solid and the liquid 
near melting differ by a few percents only. Developing 
such a theory requires a better qualitative understand- 
ing of the metastable phases of the model. It is clear 
that the overheated solid becomes unstable at some finite 
temperature. It is not clear however whether the over- 
cooled liquid becomes unstable at some finite tempera- 
ture (like water) or exists all the way down to T = as a 
metastable state. Despite its limited precision, the gaus- 
sian (Hartrce - Fock) variational calculation, is usually a 
very good guide to the qualitative features of the phase 
diagram. While such a calculation in the liquid was per- 
formed quite some time ago [18], on the solid side a more 
complicated one sampling inhomogeneous states was per- 
formed recently [9]. The results are as follows. The solid 
state is the stable one below the melting temperature, be- 
comes metastable at somewhat higher temperatures and 
is destabilized at cit ~ —5. The liquid state becomes 
metastable below the melting temperature, however in 
contrast to the solid, it does not loose metastability all 
the way down to T = and the excitation energy ap- 
proaches zero. 

Meanwhile, similar qualitative results have been ob- 



tained in a different field of physics. A variety of analyt- 
ical and numerical methods [22] have indicated that liq- 
uid (gas) phase of the classical one component Coulomb 
plasma also exists as a metastable state down to T = 
with energy gradually approaching that of the Madclung 
solid and the excitation energy diminishing. We specu- 
late that the same phenomenon would appear to happen 
in any system of particles interacting via long range re- 
pulsive forces. In fact the vortices in the London approx- 
imation resemble repelling particles with the force even 
more long range than the Coulomb. In light of this impe- 
tus to consider the above scenario in the vortex matter, 
we provide both theoretical and phenomenological evi- 
dence that the above picture is a valid one. 

Assuming the absence of singularities on the liquid 
branch allows to develop a sufficiently precise theory of 
the LLL GL model in a vortex liquid (even including an 
overcooled one) using the Borel - Pade (BP) [23] method 
at any temperature. After clarifying several issues that 
prevented its use and acceptance previously we then com- 
bine it with the recently developed LLL theory of solids 
[21] to calculate the melting line and the magnetization 
and specific heat jumps along it. Early on, Ruggeri and 
Thouless [18] tried unsuccessfully to calculate the spe- 
cific heat by using BP because their series was too short. 
Subsequent attempts to calculate the melting point by 
using BP also ran into problems. Hikami, Fujita and 
Larkin (18) attempted to find the melting point by com- 
paring the BP (liquid) energy with the one loop solid 
energy and, in doing so, obtained the melting temper- 
ature a™ = —7. However, their one loop solid energy 
was incorrect and, in any case, the two loop correction 
is necessary. As demonstrated in the following, the BP 
energy combined with the correct two loop solid energy 
computed recently not only gives a™ = —9.5, but also 
allows to obtain a wide range of quantitative predictions 
within the model. 

Now we present the solution of the LLL GL model. 
The liquid LLL (scaled) effective free energy (of the 
scaled model defined in eq.(2)) is written as fu q — 
4e 1 / 2 [l + g(x)] . The function g can be expanded as 
9 ( x ) = c nX n , where the high temperature small pa- 
rameter x = ^£~ 3 ^ 2 is defined as a solution of the gaus- 
sian gap equation, e 3 / 2 — aye 1 / 2 — 4 = for the excitation 
energy . The coefficients can be found in ref. [24,25] . We 
will denote gu{x) by the [k, k — 1] BP transform of g(x) 
(other BP approximants clearly violate the correct low 
temperature asymptotics) . The BP transform is defined 
as Jq 00 g' k (xt) exp (— t) dt where g' k is the [k, k — 1] Pade 
transform of Y^=i For k — A the liquid energy 

already converges. We used in this work k = 5 to achieve 
the required precision (<~ 0.1%). The liquid energy com- 
pletely agrees with the optimized gaussian expansion re- 
sults [9] above its radius of convergence at ax = — 5. In 
addition similar results we obtained in the 2D GL model 
agree with existing Monte Carlo simulations [15]. The 
solid effective energy to two loops is [21]: 
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Comparing the solid energy to that of the liquid (in- 
set in Fig.l), reveals that a™ = —9.5. This is in ac- 
cord with experimental results. As an example, on Fig. 2 
we present the fitting of the melting line of fully ox- 
idized YBa 2 Cu 3 O r [7] which gives T c = 88.2, H c2 = 
175.9, Gi = 7.0 1(T 5 ,k = 50. Melting lines of optimally 
doped untwinned [3] Y Ba 2 Cu^0 7 -s and DyBa 2 Cu^O'j 
[26] are also fitted extremely well. For example, for 
Y BaiGuzOis in [3], the melting line fitting gives T c = 
93.07, H c2 = 167.53, = 1.9 10~ 4 ,k = 48.5 (see also 
rcf. [27] ). 

3D Monte Carlo simulations [14] are not precise enough 
to provide an accurate melting point since the LLL 
scaling is violated. One gets different values of a™ = 
— 14.5,-13.2,-10.9 at magnetic fields 1,2,5T respec- 
tively. This violation is perhaps due to a small sample 
size (of order 100 vortices). The situation in 2D is bet- 
ter since the sample size is much larger. We performed 
similar calculation in the 2D LLL GL model and found 
that the melting point is a™ = —13.2 . It is in good 
agreement with the MC simulations [15]. Phenomcnolog- 
ically the melting line can be located using Lindcmann 
criterion or its more refined version using Debye - Waller 
factor [28]. The more refined criterion is required in the 
case of YBCO since vortices are not point-like. Numer- 
ical investigation of the Yukawa gas [28] indicated that 
the Debye - Waller factor e~ 2W (a ratio of the structure 
function at the second Bragg peak at melting to its value 
at T = 0) is about 60%. We get using methods of ref. 
[29] e- 2W = 0.59 for a?? = -9.5. 

The scaled magnetization is defined by m (ax) = 
~d^tpf e ff At the melting point a™ = —9.5 the 

magnetization jump ratio defined by AM divided by the 
magnetization at the melting on the solid side is found 
to be equal to 
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This prediction is compared on Fig. 2 (the upper in- 
set) with the experimental results on fully oxidized 
YBa 2 Cu 3 0i (ref. [7], rhombs) and optimally doped 
Y Ba 2 CuzO'j-s (rcf. [5], stars) . These samples proba- 
bly have the lowest amount of disorder which is not in- 
cluded in the calculations. From the model, we caculate 
the specific heat jump ratio at the melting 



In addition to describing the melting, we present here 
an example of quantitative results that are obtained us- 
ing the present approach - the magnetization curves. Our 
LLL magnetization curve coincides with the LLL Monte 
Carlo result of rcf. [14] (which is very accurate since the 
LLL scaling is obeyed) to the precision of MC. How- 
ever in experiments away from the melting line higher 
Landau levels (HLL) are no longer negligible. Naively 
in vortex solid when the distance from the mean field 
transition line is smaller than the inter - Landau level 
gap, 1 — t — b < 26, one expects that the higher Lan- 
dau modes can be neglected. More carefully examining 
the mean field solution reveals that a weaker condition 
1 — t — b < 12b should be used for a validity test of the 
LLL approximation [30] in vortex solids. In vortex liq- 
uid one has to go beyond the mean field to estimate the 
HLL contribution [31]. In 3D Lawrie in [31] calculated 
the excitation energy in the framework of the gaussian 
(Hartrec - Fock) approximation. The excitation gap is 10 
times smaller than inter Landau gap for fields in a wide 
range around melting line for fields larger than 0.1T in 
YBCO. Therefore in the range of values of the inter- 
est in the present paper the LLL contribution should be 
dominant. Experimentally it is often claimed that one 
can establish the LLL scaling for fields above 3 T (see, 
for example, ref. [32]). 

The theoretical expressions we use are the LLL contri- 
bution to the magnetization plus the corrections due to 
HLL calculated in gaussian approximation. The results 
are compared on Fig. 3 with the experimental magnetiza- 
tion curves of ref. [7]. We use the parameters from the 
fitting of the melting curve (see Fig. 2) . The agreement 
is fair at higher fields, while at low magnetic fields the 
higher Landau levels theory beyond gaussian approxima- 
tion is required. The LLL scaling has a limited validity 
away from melting line. 

To summarize, the problem of the quantitative descrip- 
tion of melting of the vortex lattice in the framework of 
the LLL Ginzburg - Landau approach is solved. The re- 
sults for melting line, magnetization jump and specific 
heat jump are in good agreement with experiments and 
MC simulations. This is the first quantitative theory of 
the first order melting of any kind to our knowledge. We 
believe that similar methods can be applied to other sys- 
tems undergoing the first order melting transition. 

We are grateful to T. Nishizaki and A. Junod for pro- 
viding details of their experiments. The work was sup- 
ported by NSC of Taiwan grant NSC#90-2112-M-009- 
039. 
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It is compared on the lower inset on Fig. 2 with the ex- 
perimental values of ref. [3] (using the fitting parameters 
given above). 
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FIG. 1. Free energy of solid (line) and liquid (dashed line) 
as function of the scaled temperature. The solid line ends at 
a point (dot) indicating the loss of metastability. Inset shows 
a tiny difference between liquid and solid near the melting 
point. 

FIG. 2. Comparison of the experimental melting line for 
fully oxidized Y Ba^CuzOi in ref. [7] with our fitting. In- 
set on the right shows the relative universal magnetiza- 
tion jump of 1.8 percent (line) and experimental results for 
fully oxidized Y BC12CU3O7 in ref. [7] (rhombs) and optimally 
doped untwinned Y BaiCuzO-j-s in ref. [5] (stars). Inset on 
the left shows the relative nonuniversal specific heat jump 
(line) and experimental results for optimally doped untwinned 
YBa 2 Cu 3 7 -s in ref. [3]. 

FIG. 3. Comparison of the theoretical magnetization 
curves (lines) of fully oxidized Y Ba2CuzOi utilizing parame- 
ters obtained by fitting the melting line on Fig. 2 with torque 
magnetometry experimental results [7] (dots). Arrows indi- 
cate melting points while at low magnetic field the experi- 
mental data start from the point in which the magnetization 
is reversible indicating low disorder. 
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